Conductivity anisotropy in the antiferromagnetic state of iron pnictides 
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Recent experiments on iron pnictides have uncovered a large in-plane resistivity anisotropy with 
a surprising result: the system conducts better in the antiferromagnetic x direction than in the 
ferromagnetic y direction. We address this problem by calculating the ratio of the Drude weight 
along the x and y directions, D x /D y , for the mean- field Q = (tt, 0) magnetic phase diagram of a 
five-band model for the undoped pnictides. We find that D x /D y ranges between 0.2 < D x /D y < 1.7 
for different interaction parameters. Large values of the orbital ordering favor an anisotropy opposite 
to the one found experimentally. On the other hand, D x /D y is strongly dependent on the topology 
and morphology of the reconstructed Fermi surface. Our results point against orbital ordering as 
the origin of the observed conductivity anisotropy, which may be ascribed to the anisotropy of the 
Fermi velocity. 

PACS numbers: 75.10.Jm, 75.10.Lp, 75.30.Ds 



Understanding the magnetic state of iron pnictides is 
probably the starting point to understand the origin of 
superconductivity in these compounds. The magnetic 
state is metallic, in contrast to the Mott-insulating be- 
havior of cuprates, and it presents Q = (jr, 0) columnar 
ordering: antiferromagnetic in the x direction and fer- 
romagnetic in the y direction. The measured magnetic 
moment is unexpectedly small and a structural transition 
occurs at or very close to the magnetic transition pQ. 

Recent experiments [2H4] have measured an anisotropic 
resistivity in the magnetic state with unexpected results: 
The resistivity is smaller in the antiferromagnetic x di- 
rection than in the ferromagnetic y direction. A priori 
the opposite would be expected. The anisotropy extends 
to finite frequencies as seen in optical conductivity ex- 
periments [5j [6] . Surprisingly, the resistivity anisotropy 
grows with doping and is largest close to the border of 
the antiferromagnetic phase when the magnetism and the 
structural distortion are weaker. The anisotropy can- 
not be explained by the small changes in the lattice con- 
stants produced by the structural transition as the lattice 
constant is larger in the direction with larger conductiv- 
ity [7j. It can neither be ascribed to the scattering rate 
of the carriers which has been observed to increase along 
the x direction and decrease along the y direction in the 
magnetic state [5] , as expected on the basis of scattering 
by spin fluctuations [8]. 

The anisotropy of the magnetic state has also been re- 
ported in other experiments: (i) Neutron scattering mea- 
surements [9 were described with anisotropic exchange 
constants; (ii) angle resolved photoemission spectroscopy 
(ARPES) finds a predominant zx orbital component at 
the Fermi surface at T [TJJ; and (iii) spectroscopic imag- 
ing scanning tunneling microscopy experiments [11 were 
interpreted in terms of one-dimensional bands and of 
electronic nanostructures aligned along the x-axis. 

The origin of the anisotropy is at the center of a hot 
debate. Orbital ordering (OO) in the d zx and d yz orbitals 



has been proposed [5| [8^ [TTMT8] and found in the magnetic 
state in several theoretical approaches [l3j[T4l[T9l[20]. On 
the other hand, the (71", 0)-antiferromagnetism breaks the 
tetragonal symmetry and could induce anisotropy with- 
out invoking other electronic sources [4j |2lJ [22] . Within 
this context it is crucial to uncover the origin of the ex- 
perimentally observed anisotropics. 

In a former work we have calculated the mean field 
phase diagram of the magnetic state [19 for a five or- 
bital model [23]. Two different metallic states have been 
found: the so-called high moment (HM) and low moment 
(LM) states. In the HM state all the orbital magnetiza- 
tions point in the same direction leading to a large lo- 
cal magnetic moment, except for the smallest values of 
U which give a magnetic solution [see Fig. flta)]. On 
the contrary, the LM state consists of antiparallel orbital 
magnetizations which add up to a small total magneti- 
zation. This state, stabilized by the anisotropy of the 
interorbital exchange interactions, has also been found 
in other mean-field [24] and LDA+U [25] calculations. 

OO defined in terms of a different occupation of d yz 
and d zx , n yz — n ZXl is present in both the LM and HM 
states [see Fig. []Jc)]. The orbital polarization of the 
Fermi surface shows larger zx orbital component around 
r [19] in agreement with experiments [10 . Whether the 
zx polarization of the Fermi surface requires strong OO, 
i.e. finite n yz — n ZXl or it is just a consequence of the or- 
bital reorganization at the Fermi surface due to the cou- 
pling between spin and orbitals is currently under discus- 
sion [H[22]. Furthermore, we estimated the anisotropy of 
the exchange interactions by mapping the multi-orbital 
mean field state to a Heisenberg model, and found a large 
anisotropy only in the LM state showing no correlation 
with the magnitude of the OO [T9] . 

In this letter we calculate the ratio of the Drude 
weight D x /D y and analyze how it correlates with the 
presence of OO for the metallic region of a five-orbital 
Q = (7r, 0) mean-field phase diagram to unveil the possi- 
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FIG. 1: (Color online) (a) U vs J/U magnetic phase diagram superposed to the magnetization showing the paramagnetic 
(black), high moment (HM) and low moment (LM) magnetic phases. The magnetic moment reaches values as high as 3.4 
\±b- The color scale emphasizes the region with magnetic moments smaller than 2 \±b- (b) Drude ratio with largest values of 
D x /Dy > 1 corresponding to small values of the magnetization, (c) Orbital ordering n yz — n zx . It correlates well with the 
difference between the partial density of states of d yz and d zx at the Fermi level, see supplementary material. 



ble source of the transport anisotropy observed in pnic- 
tides [IHS] . We find a wide range of values for this ratio 
(0.2 < D x /D y < 1.7) depending on the value of the in- 
teraction parameters. Larger 00 appears together with 
larger conductivity in the ferromagnetic direction, there- 
fore it cannot be responsible for the observed anisotropy. 
We show that the magnitude of the Drude anisotropy 
is related to the topology and morphology of the Fermi 
surface. The observed conductivity anisotropy may be 
ascribed to the Fermi velocity, as in Ref. [4]. 

The 5-orbital interacting model Hamiltonian and 
mean- field U vs J/U Q = (tt, 0) magnetic phase di- 
agram were described in detail in Ref. [19 . As de- 
tailed in Ref. [23], both direct Fe-Fe and indirect (via 
As) hoppings, calculated within the Slater-Koster frame- 
work [26 , determine the magnitude of the hopping ampli- 
tudes which depend on the angle a formed by the Fe-As 
bonds and the Fe-plane. Here we use a = 35.3° cor- 
responding to a regular As tetrahedra. Energy units 
are ~ leV. The interactions include the intraorbital 
Hubbard [/, the Hund's coupling J, and the interobital 
U' = U — 2J. The Hamiltonian is solved at the mean- field 
level keeping only the spin and orbital-diagonal average 
terms. In Fig. [IJa) we reproduce the phase diagram in 
the metallic region (U < 2.8 eV) superposed to the re- 
sulting magnetization. As discussed above, HM and LM 
solutions arise for high and low values of J and both show 
OO, see Fig. flic). OO observed here has the same sign as 
the one obtained in LDA+U [13^ and LDA+DMFT [16]. 

The Drude weight involves the carrier velocity, the 
spectral weight, and the scattering rate. Experimentally, 
the scattering rate is anisotropic, but with an anisotropy 
opposite to that of the resistivity [5]. Here we neglect 
the influence of the scattering rate. This approximation 
allows us to determine the role of the Fermi velocity in 



the observed transport anisotropy. We calculate 
Z Kn v 2 x (k,n)g(k,n)S(e n (k)-E F )) 



D x /D y 



E k ,« w y( k > n )5(k, n)<5(e„(k) - E F ) 



(1) 



with Vi the velocity in the direction i = x,y, n the band 
index, g the spectral weight, e n the energy of the band n, 
S the Dirac function, Ep the Fermi energy, and the sum 
in k restricted to the Fermi surface. The result for the 
Drude weight ratio is shown in Fig. mb). In the param- 
agnetic (PM) region, black region in the magnetic phase 
diagram Fig. flla), D x /D y = 1 as expected (discrepan- 
cies along the boundary region might be due to numer- 
ical precision). The magnetic phases show anisotropy 
in the Drude weight with ratios 0.2 < D x /D y < 1.7. 
D x /D y > 1, consistent with the experimentally observed 
transport anisotropy [2HI] , is observed within a great part 
of the LM region for intermediate values of the on-site in- 
teraction well suited for describing pnictides. In the HM 
phase, D x /D y > 1 is restricted to a region close to the 
PM phase, corresponding to the smallest values of the 
magnetic moment. Interestingly, in a wide part of the 
region with D x /D y > 1 an increase of the anisotropy 
correlates with a decrease of the magnetic moment. On 
the other hand D x /D y > 1 happens to be when the OO 
is small (\n yz — n zx \ < 0.05). Note that although there 
is not a perfect one to one correspondence between the 
anisotropy of the Drude weight [Fig. [IJb)] and the OO 
[Fig. flic)], D x /D y decreases with increasing OO in most 
part of the phase diagram. That is, OO seems to favor 
a transport anisotropy opposite to the one reported ex- 
perimentally [4]. This finding rules out OO as the origin 
of the observed resistivity anisotropy. 

We now analyze the Drude weight results on the light 
of the shape of the reconstructed Fermi surface and the 
velocity vector on it [4 . These are represented in Fig. [2] 
in the extended (Fe) Brillouin zone for different points 
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FIG. 2: (Color online) Fermi surfaces (left) and velocities 
(right) at different points of the phase diagram in the ex- 
tended (Fe) Brillouin zone. The Fermi surfaces are weighted 
by their spectral weight. The corresponding U and J val- 
ues, as well as the obtained anisotropy of the Drude ratio, 
magnetic moment in Bohr magnetons and the magnetic state 
appear on top of each figure. 



of the phase diagram corresponding to different values of 
the ratio D x /D y . We have chosen Fermi surfaces repre- 
sentative of a wide area of the phase diagram, while other 
topologies have also been found. We observe a strong de- 
pendence of the morphology and topology of the Fermi 
surfaces on the interaction parameters. These variations, 
which have a correspondence with the calculated Drude 



weight ratio, can be understood starting from the PM 
Fermi surface in Fig. [3la) and the plot of the bands in 
Fig. [3Fb). The bands are shown for the PM state (in 
black) and the LM state (in red) for U = 2.2 eV and 
J/U = 0.07. 

The PM Fermi surface has hole pockets in T and M, 
and electron pockets in X and Y. In the magnetic states, 
a gap opens at M and Y near Ep but, for small values of 
the magnetic moment, there are still pockets in M —> Y 
though considerably reduced in size with respect to the 
PM phase, see Fig. [3] This is illustrated in Fig. [2|a) for 
U = 2.2 eV and J/U = 0.07 (corresponding to the LM 
state) and is also representative of the HM state with a 
small value of the magnetization (region close to the PM 
phase). This topology of the Fermi surface corresponds 
to the largest values of the anisotropy D x /D y > 1.2. A 
visual inspection of the velocity vectors at kp in Fig.^b) 
shows that the v x component of the velocity dominates 
in this case, with a large contribution from the pockets at 
M — »> Y. D x /D y > 1 appears also at points of the phase 
diagram with the Fermi surface represented in Fig. ^1 (c) 
for U = 1.9 eV and J/U = 0.15. The main difference 
between Figs. [2|a) and (c) is the disappearance of the 
pockets at M — >• Y in the latter case while the shape of 
the pockets at F and along F — >• X is maintained. The 
electron pockets along F — > X are associated to two Dirac 
cones just below Ep as shown in Fig. [3j The morphol- 
ogy of the Fermi pockets at F —t X in Figs. [2|a) and (c) 
also supports D x /D y > 1 as evidenced in Fig. J2|d). It 
is not straightforward to make a connection between the 
anisotropy of the Drude weight and the orbital content 
of the Fermi pockets. The pockets close to Y show zx 
and xy polarization [23] , while in the Dirac cone electron 
pocket the v x contribution comes mostly from the d yz po- 
larized segment closer to F [19]. In fact, the anisotropy 
of the Drude conductivity cannot be ascribed to the or- 
bital reorganization at Ep, whose dependence on U and 
J follows that of the OO, see supplementary material. 

Now we turn to Fig. ^e) (for U = 2.6 eV and J/U = 
0.05). The topology of the Fermi surface is maintained 
with respect to Fig. p[c) but the morphology is different: 
the pockets are more isotropic and those close to X have a 
slightly stronger v y component rendering D x /D y < 1 [see 
Fig.^f)]. Finally, Fig. ^g) (for U = 2.8 eV and J/U = 
0.145) corresponds to a strong anisotropy D x /D y <C 1, 
opposite to the one observed experimentally. The pockets 
along r — >> X have disappeared. Two ellipsoidal pockets 
with a strong v y component close to this direction are 
observed instead, see Fig. [2^f ). 

The Fermi surfaces in Fig. [2|a) and (c), once folded 
to the FeAs Brillouin zone and taken into account the 
existence of magnetic domains, resemble the flower-like 
Fermi surfaces observed in ARPES in the magnetic 
state [27[ [28]. From the available experimental results 
we cannot conclude on the presence or absence of Fermi 
pockets close to Y. 
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FIG. 3: (Color online) (a) Fermi surface in the paramagnetic state for U,J = 0. (b) Band structure for U = 2 eV and 
J/C/ = 0.07 corresponding to the LM state, in red. The linewidth gives a measure of the spectral weight. The paramagnetic 
bands, in black, are shown as a reference. The position of the Dirac cones along YX is marked with arrows. 



In conclusion, we have found that the anisotropy of 
the Fermi velocity which is strongly related to the mor- 
phology and topology of the reconstructed Fermi surface 
results in large variations in the ratio of the Drude weight. 
Depending on the shape of the Fermi surface the resis- 
tivity can be larger in the antiferromagnetic x or the 
ferromagnetic y direction. Larger orbital ordering cor- 
relates with a resistivity anisotropy opposite to the one 
observed, i.e. larger conductivity in the ferromagnetic 
y direction. A larger conductivity in the antiferromag- 
netic x direction, consistent with the experimentally ob- 
served conductance anisotropy, corresponds to Fermi sur- 
faces whose shape resembles those observed experimen- 
tally. Values of the anisotropy D x /D y > 1.2 are achieved 
when the Fermi surface has pockets at Y in the extended 
Fe Brillouin zone. Pockets around Y are also predicted 
in ab-initio [21] and mean field calculations [29l [30] . The 
electron pockets originating from the Dirac cones can also 
govern the anisotropy when D x /D y > 1. The region of 
the phase diagram with D x /D y > 1 shows a small mag- 
netic moment, as observed experimentally. 
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